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Abstract. The non infinite divisibility property is proved for a class of photon counting 
probability distributions attached to Euclidean Landau levels and a decomposition of 
their corresponding laws is explicitly obtained. 



1. Introduction 

It has been known that the photoelectric process is essentiahy quantum in nature [T]. Its 
proper explanation involves the particle properties (photons) of the electromagnetic field 
and is thereby able to explain the statistical photoemission. However form a semiclassical 
point of view the incidence of the electromagnetic field may be regarded as modifying the 
statistics of the counter. A quantitative relationship between intensity and the probability 
of photoemission, and hence counts, may be deduced from the conventional time-dependent 
perturbation theory of quantum mechanics. Indeed, if we suppose temporarily that there 
are no random fluctuations in the intensity of the electromagnetic field and we assume that 
the probability of counts in distinct time intervals are statistically independent and that 
the average number of photons per unit time is a constant A then the resulting appropriate 
probability distribution X is the familiar Poisson law with parameter A (|21 pp. 16-19]). 

This probability distribution arise also form the standard coherent states of the har- 
monic oscillator in the Fock-Bargmann representation involving the Bargmann space of 
entire Gaussian square integrable functions (O p. 70- 72]). The latter one can be realized 
as the eigenspace of the Hamiltonian operator with uniform magnetic field normal to the 
plane, which is associated with the lowest Landau level, say, eo Moreover, to any 
eigenspace of this Hamiltonian, associated with an mth Landau level, say €m, correspond 
a class of generalized coherent sates giving arise to random variable X counting the number 
J = 0, 1, . . . , of photons with probabilities [3]: 



:i.i) 



where mAj := min {m,j), mVj := max {m,j) and (.) denotes the Laguerre polynomial 
[B]. The expression p.l|) may also appear when calculating the probability transitions in 
some special cases of molecular vibrations, which are called Pranck-Condon factors [7]. 
From p.ip . it is clear that X reduces to Poisson law V (A) for m = 0, which is infinitely 
divisible [8j in the sense that for any integer n it decomposes as a sum of independent 
identically distributed random variables Xi as X = Xi + - ■ ■ + Xn- Therefore, one naturally 
wonders if X is still infinitely divisible for general m = 1, 2, 3, • • • . 

In this, paper we shall prove that the infinite divisibility breaks down when m > 1. More 
precisely, if we denote by x^^\ 1 < k < m, the ordered zeros of the Laguerre polynomial 

(x), then X is not infinitely divisible whenever < A < \x^i"^ or jx^™'^ 
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integer k. We shall also give a decomposition of the law of X as 

li~ = V (A) *5m*(Y. \lk (A, m) {5k + , (1.2) 

where 5a are the Dirac mass and the coefficient 7^ (A, m) are given explicitly as polynomials 
functions of degree m with respect to A. 

The paper is organized as follows. In Section 2, we review briefly the coherent states 
formalism used to construct photon counting probability distributions. While in Section 3, 
we recall some needed facts on the so-called generalized Bargmann spaces on the complex 
plane. Section 4 is devoted to recover the Poisson probability distribution as photon 
counting distribution associated with the Bargmann space and to discuss related known 
properties. In Section 5, the construction of the generalized probability distribution X 
based on a class of generalized coherent are recalled and relevant characteristic properties 
are also presented. In Section 6, we prove the non infinite divisibility of X and we establish 
a decomposition of its probability measure. 



2. Coherent states 

Here, we review a coherent states formalism starting from a measure space A" "as a set of 
data" as presented in [3]. Precisely, let {X, a) be a measure space and let C L'^{X, a) be 
a closed subspace of infinite dimension. We now suppose that £ is a self-adjoint operator in 
a separable Hilbert space, say, "H, of quantum states. Let us assume that the spectrum of 
£, has only a discrete component. Normalized eigenstates of C are denoted by {| cpj >}^o 
and they form an orthonormalized basis of H. Next, suppose that this basis is in one-to- 
one correspondence with an orthonogonal basis {^j}J^Q of A'^. Furthermore, and this is a 
decisive step in the wave-packet construction, we assume that 

00 

0<Mix) = J2pf\'^,{x)\^<+oo, (2.1) 

j=0 
1 1 2 

almost everywhere on X, where pj := \\^j\\i2(^;Y ^^y Then, the states in Ti 

|x>:=(AA(x))-2^^|<^,> (2.2) 
i=o VPj 

are coherent states labeled by points x €z X. By this definition, the nornmalization 
< X \ X >= 1 is ensured because of the orthonormahty of the set {| >} and the 

presence of the normalization factor (M {x))~^ . The resolution of the unity !■}{ in T-L holds 
by virtue of the orthonormahty of the set Indeed, we have, in terms of Dirac's 

bra-ket notations: 

/ \x><x\M{x)a {dx) = ^\(l)j X (t^k I / {x) {x)o- (dx) (2.3) 
Jx .^^ Jx 

= 5] I 0,- 1=1^, (2.4) 

j 

where we have used the orthonormahty of the <I>^s and the resolution of the unity in % 
obeyed by the orthonormalized basis {| >}. A statistical aspect of these coherent states 
is that we have, for almost each x, a discrete distribution 
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VP] 
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AB = -^,r7^ + Bz]—. (3.4) 



which could be considered as concerning experiments performed on the system to measure 
spectral values of the "quantum observable" 

3. The generalized Bargmann spages (C) 

In suitable units and up to additive constant, the motion of a charged particle in a 
constant uniform magnetic field in is described by the Schrodinger operator 

Hb = -\ ((9., + ^By,)" + (dy^ - ^Bx,)"^ - \ (3.1) 

acting on L^(M^, d/x), d[i being the Lebesgue measure on R^, i3 > is the strength of the 
magnetic field. We identify the Euclidean space with C in the usual way to use the 
complex structure. The operator Hb in (|3.1|) can be represented by the operator 

AB:=ei^N'i/Be-i^l^l'. (3.2) 

According to ()3.2|) . an arbitrary state (\) of L^(M^, d/i) is represented by the function Q [0] 

of L2(Ce--^l^l'd^) defined by 

Q[^](z) :=e5-S'^'V(2), ^eC. (3.3) 

The unitary map Q is called ground state transformation. Explicit expression for the 
operator in ()3.2p is: 

From now on, we take i? = 1 in ()3.4p and we consider the magnetic Laplacian A := Ai 
acting in the Hilbert space L^(Ce~l^l d^) of complex-valued Gaussian square integrable 
functions. The spectrum of the operator A consists of eigenvalues of infinite multiplicity 
of the form (Landau levels): 

€m-=m, m = 0, 1,2, • • • . (3.5) 

The corresponding eigenspaces 

Al (C) := {99 G l2 (c, e-l"l'd^) , Aif = em^} (3.6) 

are called generalized Bargmann spaces on C. For each fixed m, the space in (j3.6p admits 
an orthogonal basis 

$„,(z) = (-1)™^^- (m A j)! kl'""^'' e-('"-^)-8^4';:7^"l) (iz^) (3.7) 

with the norm square given by 

Pm,j ■= ll^m.jll^ = 7rm!j!. (3.8) 
Finally, let us note that the eigenspaces in (|3.6p satisfy 

(C, e-l^l'd^) = AliC) ® AliC) e ■ ■ ■ e AliC) e ■ ■ ■ . (3.9) 

By [1], it was proved that the eigenspace Aq (C) corresponding to the lowest Landau level 
eg = coincides with the well known Bargmann space (C) consisting of entire functions 

in L2 (c, e-l^l'd^) . That is, A^ (C) = ^ (C) . 
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4. The Poisson law V (A) and the eigenspace Aq (C) 

In this section, we recover the definition of the Poisson distribution as a photon-counting 
probability distribution attached to the lowest Landau level eg. We will also recall some 
known properties of this law. 

In order to apply the above coherent states formalism, we first list the following elements. 

• {X,a) = ^C,e~l^l with a := e"'^' d^, /i is the Lebesgue measure on C. 

• := Al (C) C L2 (C,e-l^l'd/i) is the eigenspace of A corresponding to the 
lowest Landau level eo or the Bargmann space F (C) . 

• Mq{z) = 7r~"^el^l is a normalization factor 

• An orthonormal basis of (C) is given by 

$o,i(^):=^, i = 0,1,2,... (4.1) 

• The Hilbert space T-L can be taken as L? (M) on which the quantum harmonic 
oscillator C = —d'^ + is acting. 

• are elements of an orthonormal basis of H consisting of eigenstates of the 
operator C. 

hoo 



>:= (AAo(z))-^ ^ ^ I >, zGC. (4.2) 



In view of (j4.1|) the probability of finding j photons in the coherent state | 0^ > is given 
according to (|2.5|) by the squared modulus of {(j)z-, 4>j) as 

Pr (X = j) = p, (|z|^) = \{(t>zA,)H? = ^! (^-^^ 

which is recognized as a standard Poisson probability distribution with a mean value 
A := \z\ . The parameter A represents the average number of photons in the coherent 
states \ (pz >■ 

We now recall some known properties of X ^ V (A). 

1. Let [x] denotes the greatest integer not exceeding x. Then cumulative function 
X r^V{X) reads 

(^) = E P. (A) = E ^ = rjryr (1 + [x];X) , (4.4) 
i=o i=o ^^J- 

where T {1 + q;u) denotes the incomplete gamma function [6]. 

2. The characteristic function (c.f) of the random X V (A) is defined by the formula 

+ 00 

^xiu) = E (e»^) = e^^yj (A) , (4.5) 

where m is a real number, i = is the imaginary unit. Explicitly, 

(^i(7x) = exp(A(e™-l)). (4.6) 

3. The mean and the variance are given by 

E{X) = Var{X) = X. (4.7) 
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4. Note also that (14.61) can be also written as 



V>°xiu) = (^^x^{u)^ , (4.. 



where 



"fx 1 [u) = exp 



^(e»-l)) (4.9) 



is the characteristic function of a Poisson probability distribution with parameter X/n. 
This means the well known infinite divisibility of V (A) . 

5. The generalized Poisson law V (A, m) and the eigenspace A^^ (C) 

We review briefly the construction of coherent states attached to Landau levels as [S]. 
In order to adopt (j2.3|) . we first list the following elements. 

• {X,a) = (C,e-l^l'(i/i) 

• A"^ := A"^ (C) C ^C,e~l^l dfij is the generalized Bargmann space in (j3.6|) . 

^ are the eigenfunctions of A given in p.7p . 

• Pm,j = vrm!j! as in (j3.8|) . 

• M{z) = Tr'^e'^l is a normalization factor. 

• {(pj}^^ are elements of an orthonormal basis of a Hilbert space such that dim 
n = dimAl, (C) = CO. 

For each fixed m = 0, 1, 2, • • • , we define a set of coherent states according to (|2.3|) by 

+00 I / \ 

I ^.,m >= (AA(Z))-^ Y: I > • (5-1) 

The system being prepared in the coherent state ()5.ip . then the probability of finding j 
photons is given according to (12. 5p by the squared modulus of the scalar product 

(</>.,n^,<Ai>« = (AA(z))-^^!^. (5.2) 

This probability is of the form 

for varying j = 0, 1, 2, . . . ,with A = In this case, the random variable X counting 
the number j of photons was denoted X ~ "P (A, m) and called in [T] a generalized Poisson 
distribution attached to space A'^ (C) or to the mth Landau level. We now will summarize 
some of its principal characteristic quantities as follows. 

1. We start by the distribution (cumulative) function (d.f): 

[x] 

Kix) =J2pj{X,m) , xG[0,+oo), (5.4) 
i=o 

where [x] denotes the greatest integer not exceeding x. We can make use of the Charlier 
polynomials ([6J): 

Q,il,X)=qlL^-HX), (5.5) 
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where q is the degree of Laguerre polynomial Lq '^^ (A) and (/ — q) the order to rewrite 
(1^ as 

-A H 



F^{x) = —Y^——{QmM{{myj)-\))\ xeK. (5.6) 
m! V. 

2. The characteristic function (c.f) of the random X ^ V {\, m) is defined by the 
formula 

/ ~\ 

VP^(n) = E e-^ = ^ e*^>,- (A, m) , (5.7) 
^ ^ i=o 

where u is a real number, i = is the imaginary unit. Explicitly, we 

^|(n) = exp(A(e*"-l))e^'""L(^)(2A(l-cosn)), ue^. (5.8) 

3. It is known that if A'I~{u) := (p~(—iu) exists, at least in a neighborhood of n = 0, 
then M~{u) is called the moment generating function (m.g.f) of X. In this case all moments 
of X exists and can be obtained by higher differentiation of M~(m) as 

In our case, the moment generating function is of the form 

M^(ti) = exp (A (e^* - 1) - mu) Lm (2A(1 - cosh-u)) . (5.10) 

4. The computation of the first and second moment gives the mean and the variance 
of X as follows 

e(x) = A + m,yar (x) = A(l + 2m). (5.11) 

Remark 5.1. For m = 0, the expression in (jS.lip reduces to (j4.7p the which are charac- 
teristic of V (A) . 

6. Analysis of X ~ P (A, m) 

Recall that a random variable Y is infinitely divisible if for any integer n it decomposes 
as a sum of independent identically distributed (iid) random variables Y^j^^, 1 < k < n 
(called infinitesimal parts, see [8]): 



y = T.y!/l (6-1) 



■ l/n- 
k=l 

In particular, the characteristic function does not vanish and this is only a necessary 
condition for infinite divisibility on the real line; (see [S] for the general theory). As 
matter of fact, in our setting, for any m > 1, the parameter A must satisfy 

2A(1 - cosn) 7^ l<k<m, (6.2) 

for all ti G R, where x^™'* ,1 < k < m are the ordered zeros of Laguerre polynomial 
(x) , which are known to be simple and positive [S]. But since 

< 2A(1 - cosu) < 4A, (6.3) 

then one has to do the restriction A < jx^^^ to avoid the vanishing of the characteristic 
function when looking if X ~ P (A, m) is infinitely divisible or not. In fact, taking into 
account this discussion on the location of the parameter A with respect to zeros of Laguerre 
polynomials, we establish the following result 
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Theorem 6.1. Let x^j^\ 1 < k < m, denote the ordered zeros of the Laguerre polynomial 

Lm (x). Then, if < A, for some k, or < X < \x^^\ then X ^ V (A,m) is not 

infinitely divisible. 

Proof. If for some fc, 1 < /c < m, we have that \x^^^ < A then the characteristic function 
^xiu) in (|5.8|) will vanish at the point 

Uk ■ 



I 1 ^ C™) 
ar cos I 1 — TTT-^^fc 
2A 



(6.4) 



and then X can not be infinitely divisible. But if A < then 2A (1 — cos u) ^ x"^'' , 1 < 

k < m, for all u G M, where the ordered zeros x^^\ 1 < k < m of Laguerre polynomial 
Lm (x) are known to be simple and positive, so that the characteristic function can not 
vanish. Now, since X is discrete, we shall use for instance the following criterion for infinite 
divisibility (O p. 51]): A discrete probability distribution {pk)k£N with po > is infinitely 
divisible if and only if the logarithmic derivative of its generating function 

P{z):=Y,Pkz\ ze[0,l] (6.5) 
fc>0 

satisfies 



P(z) 



fe>0 



(6.6) 



for z S [0, 1) . In order to apply this criterion, we note first that the characteristic function 
$~ extends to an analytic function on the upper-half plane C"*" since X is valued in N. 

As a matter of fact, the exponential moment generating function of X is well-defined for 
u>0: 



E e 



exp (A (e " — 1) — mu) Lm (2A (1 — cosh n)) . 
Setting z = e~" then z E [0, 1] and the probability generating function of X reads 



P{z) = E[z 



X 



exp 



(A (z-1))z™lW _ (!_,) 



A 



Hence, a direct calculation give successively 



(6.7) 
(6.8) 



P{z) 



A + - + ^LogL^^^ 

z dz 

J m 

X + — + — J2Log 



z 



dz 



k=l 



+ -(!-.)- 



A + 



A(l 



t^izxt^^ + Kl-zY^ 
A (1 - z^) 



k=l 



ZX 



(m) 



+ x{i-zy 



X 



(m) 



2A(l-z2) 



Szxf) + A(l-z)2' 
This defines an entire function whose constant term is computed as 



A 



2 III' 

i=l 



2mA + m' 



(6.9) 
(6.10) 

(6.11) 

(6.12) 

(6.13) 

(6.14) 
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where the last equality follows from the expression of the Laguerre polynomial 

Lm{x) = f: . (6.15) 

Thus the constant term is obviously positive and we shall prove that the first derivative 
at z = is negative. By direct calculations, this derivative equals 

_ J_ |^2m A2 - 4Am2 + ^ (xj™) ) . (6.16) 

Next, we use the fact that 

^2(4""^) =m4-2^xi'")a;;") =m2(2m-l). (6.17) 

1=1 k<l 

Thus, we need to check that the polynomial 

i?(A) := 2mA^ -4Am2 + m2(2m- 1) (6.18) 

has a strictly positive sign for < A < ^x^-^^K The zeros of R (A) are 

A_ := m — -\/2rn, A+ := m + -\/2m. (6.19) 

So i? (A) > for A ^ [A_, A+] . We need to check that 

ix^™) < A_. (6.20) 

This is why we need an estimate on the "first" zero Xi^\ By ([9j, p. 125]) the following 
estimates holds for the least zero x^™^ for the Laguerre polynomial Lm^ (x) : 

(5i/2)^ (m) (a + l)(a + 3) , , 
<xr' < \ \ , m = 0,l,2,---, (6.21) 



m + (a + 1) /2 ^ ~ 2m + a + 1 ' 

where 5i stands for the least positive zero of the Bessel's function (x) . Applying this 
in the case a = 0, we obtain that 

l(™)<^/i^^ m = 0,l,2,..., (6.22) 
4 2m + 1 



Now, we can check that 



3/4 I , , , , 

' <m--V2m = A_ (6.23) 



2m + 1 2 

for m > 1. This ends the proof. □ 

Remark 6.2. Infinite divisibility behavior of the ground state probability distribution for 
some selected quantum systems in the Euclidean quantum mechanics setting was discussed 
in \W\. 



Now, to establish a decomposition for the measure of X ~ (A,m), we first need the 
prove the following result. 

Lemma 6.3. For m = 0, 1,2,---, A > Oand n G M, the Laguerre polynomial in (|5.8|) 
decomposes as 

Lm (2A (1 - cos u)) = ^° ^ ' ^' +^7k (A, m) cos ku, (6.24) 

^ k=i 
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where the coefficients are given by 



-m) 



7.(A.»):=2(-1) — , (6,25) 



Proof. We start by writing formally (6.24) then the coefficients 7^ (A, m) can be obtained 
as a Fourier transform 

1 

7fc(A,m) = — / (cos A;u) Lm (2A (1 — cos ti)) du. (6.26) 
vr 7o 

We set u = 2x, then (|6.26|) takes the form 

7fe(A,m) = — / (cos 2/cx) Lm (4A sin^ X ) (ix. (6.27) 
vr 7o ^ ^ 

Next, we make use of the formula ([11, p.225]): 

/ cosaxLl^Mfesin^xMx = -sin cos (6.28) 

JO ^ ^ a 2 2 m! 

XaFa (y-m, ^, 1, a + 1, 1 - /c, 1 + A;; 6^ 

for a = 2k, b = 4A, n = 1 and a = 0. In this particular case, the hypergeometric function 
3F3 reduces to an as 



3 



F3 [-m, ^, 1, 1, 1 - fc, 1 + k; 4A^ =2 r2 (^-m, ^,l-k,l + k; 4A^ . (6.29) 



The latter also reads 



r(l-A;)r(l + A;)(-m)^.(i) 



/ 1 \ ™ i ft-M IJ^-^^'M-"'-;,- I 2; . MAV 



. >. + fc + i)r(i-k + j) j! 

where the singular part coming form T (1 — k) has to be understood as Euler's reflection 
formula [B]: 

TT "r(l-A:)r(A:)- ^^'"^ ^ 

Therefore, 



(^), (4A)^ 



7fe (A, m) := 2 (-1)^ m ^ 1. ^ Wl ' " ^^'^^^ 
^^r(l + /c + j)r(l - A: + j) j! 

The required formula for 7^ (A, m) follows from the fact that z 1— )■ 1(2;) is meromorphic 

with poles the negative integers. □ 

Proposition 6.4. The law of the photon counting random variable X ~ P (A, m) can be 

decomposed as 

f,~ = V{\)*6^* (f;i7fc(A,m)(5fc + 5_fe) ) , (6.33) 



\k=0 ^ 



where 7^ (A,m) are gifwen by (|6.32|) . 
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Proof. The characteristic function in (|5.8|) can be sphted to a product of three pieces: 
exp (A (e*" — l)), e*™'* and Lm (2 A (1 — costi)) corresponding to three measures, namely 



exp 



(A(e™-l)) ^ P(A) (6.34) 



m 
m 



(2A (1 - cosn)) < — > ^ (A, m) (4 + <5-fc) , 

fc=o 

where P (A) is the standard Poisson distribution, 5m is the Dirac mass at the point m and 
the coefficients 7^ (A,m) are given in Lemma |6.3[ □ 

Remark 6.5. The coefficients 7^ (A, m) can take negative values. For instance when 
m = 1, the Laguerre polynomial reads 

(2A (1 - cos ti)) = 1 - 2A (1 - cos n) = (1 - 2A) + 2A cos u. (6.35) 

Here 70 (A, 1) = 1 — 2A and 71 (A, 1) = 2A. It is clear that the sign of 70 (A, 1) changes 
according to X < or X > ^. We should look very closely to sign of the coefficient 7^ (A, m) 
or to the polynomial 

n (A) :=2 r2 [-m, ^,l-k,l + k; 4A^ (6.36) 

or at least 70 (A, m) may be it could encode information about super/sub Poissonian statis- 
tics of the photon counting distribution X ~ "P (A, m) since this character depends on the 
wether A < ^ or A > |; see Proposition ?? in [5]. The classical character (Poissonian) of 
the distribution is encoded in the Poisson distribution V (A) part of the decomposition. 
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